Localization of preformed Cooper-pairs in disordered superconductors 
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The most profound effect of disorder on electronic systems is the localization of the 
electrons transforming an otherwise metallic system into an insulator. If the metal 
is also a superconductor then, at low temperatures, disorder can induce a dramatic 
transition from a superconducting into an insulating state. An outstanding question is 
whether the route to insulating behavior proceeds via the direct localization of Cooper 
pairs or, alternatively, by a two-step process in which the Cooper pairing is first de- 
stroyed followed by the standard localization of single electrons. Here we address this 
question by studying the local superconducting gap of a highly disordered, amorphous, 
superconductor by means of scanning tunneling spectroscopy. Our measurements re- 
veal that, in the vicinity of the superconductor-insulator transition, the coherence 
peaks in the one-particle density of states disappear while the superconducting gap 
remains intact indicating the presence of localized Cooper pairs. Our results provide 
the first direct evidence that the transition in our system is driven by Cooper pair 
localization. 



Although superconductivity and Anderson localization 
lead to the opposite extremes of conductivity at low tem- 
perature (T), both are due to delicate quantum effects. 
In superconductors, electrons are bound in Cooper pairs 
that condense into a zero-resistance, macroscopic, quan- 
tum state. In contrast, disorder induces quantum local- 
ization of the electron's wave-function that transforms a 
metal into an insulator with diverging resistance. 

It turns out that an increasing level of disorder can 
cause a transition from a superconductor into an in- 
sulator. Understanding how the disorder drives this 
transition is important for many quantum systems 
such as amorphous superconductors [l[ , superconducting 
nanowires [2j, high critical-temperatures superconduc- 
tors 3 and ultra-cold atomic gases Q • Furthermore, this 
transition is regarded as one of the prototypical quantum 
phase-transitions driven by disorder in a many-body sys- 
tem, a subject that acquired significant theoretical atten- 
tion recently. 0-01 • 

For a moderate level of disorder the Anderson the- 
orem [H, which is based on the Bardeen-Cooper- 
Schrieffer (BCS) theory of superconductivity [Io| . states 
that the critical temperature of superconductivity, T c , 
remains unchanged. In order to affect T c significantly a 
much stronger disorder, at a level that usually causes lo- 
calization of the electronic wave functions, is needed. As 
a result the suppression of superconductivity is theoret- 
ically expected to be accompanied by a transition to an 
insulating state [Ill4l7|. 

The route to the total destruction of superconductiv- 
ity by disorder can follow two distinct paths. In the first, 
and more obvious, path disorder-enhanced Coulomb re- 



pulsion eliminates Cooper-pairing before the onset of lo- 
calization. The ensuing poor metal becomes an insulator 
upon a further increase of disorder [f8j]. In this case one 
expects to find an insulating regime similar to that found 
in other, non-superconducting, metals driven to become 
insulators by strong disorder. 

The second, and more intriguing, path is where the su- 
perconductor itself undergoes the transition into an in- 
sulating state with only minimal sup pression of Cooper- 
pairing. For such materials, [f9r|26| there are two main 
theoretical mechanisms that have been suggested to ex- 
plain the transition [27j | . These mechanisms differ by 
their emphasis on the relative importance of processes 
on short and long scales. In non-granular materials, 
which are the focus of this paper, the short length scale is 
set by the coherence length of the superconductor. The 
first of these mechanisms attributes the superconductor- 
insulator transition (SIT) to the divergence of phase fluc- 
tuations at large scales |28[ . The predictions of this the- 
ory are in agreement with resistivity measurements in 
quench- condensed Bi films The second suggested 



mechanism [17|, |29( emphasizes the role of the fractal 
nature of the electron wave-functions at short length 
scales. Extending earlier works llHl4|. this theory pre- 
dicts that superconductivity at high disorder is main- 
tained by a fragile coherence between a small set of pre- 
formed Cooper-pairs that are characterized by an anoma- 
lously large binding energy. Consequently, in the vicinity 
of the SIT, both the insulator and the superconductor are 
composed of these preformed Cooper-pairs that either lo- 
calize, leading to an insulating state, or condense into a 
coherent zero-resistance state. 
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The existence of preformed Cooper-pairs has been pre- 
viously inferred from transport measurements findings 
such as the giant magneto-resistance peak and activated 
resistance in the ma gnet ic field-tuned SIT of Indium- 
Oxide (InO) (2(1 HHTH, or thin Titanium-Nitride 



(TiN) films [31| and, more recently, from magneto- 
resistance oscillations with half flux-quantum periodicity 
in insula ting Bi films patterned with a honeycomb array 
of holes 32 , 33| . While these observations are in agree- 
ment with the presence of preformed Cooper-pairs at the 
SIT [29[ they only constitute indirect evidence for their 
existence. 

The existence of Cooper-pairs at short length-scales is 
more directly revealed by the observation of a supercon- 
ducting gap in the density of states (DoS) [13, 0] using 
scanning tunneling microscope (STM) spectroscopy [3^ . 
In good BCS superconductors global coherence sets in 
at T c . Below T c , a DoS suppression begins below the 
gap energy, A, and the lost spectral weight appears as 
distinct peaks at A. These 'coherence' peaks are closely 
linked to the emergence of a long-range superconduct- 
ing state. In the case of preformed Cooper-pairs with- 
out global coherence, it is theoretically expected 17, 29[ 
that the spectrum will remain gapped but the coherence 
peaks will be missing. Only when global superconductiv- 
ity sets-in, the coherence peaks re-emerge. The height of 
these coherence peaks is predicted to fluctuate strongly 
from one location to another. 

In this article, we report on a systematic spectro- 
scopic study of the local properties of superconductivity 
in highly disordered amorphous InO films close to the 
SIT. Our results, obtained using an STM mounted in a 
custom-made dilution refrigerator capable of a base tem- 
perature of 50 mK, provide the first direct evidence for 
the existence of preformed Cooper-pairs along with their 
localization. 



Localization of preformed Cooper pairs 

The main feature of our results is the existence of two 
qualitatively different T-evolutions of the DoS spectra as 
illustrated in Fig. QJi and QJ). For T > T c (T c is in- 
dicated by the black dashed trace in the figures), both 
T-evolutions exhibit very similar behavior with a low- 
energy DoS suppression that deepens as T is lowered to- 
wards T c . A similar DoS suppression above T c has been 
seen in other superconductors and has been dubbed the 
pseudogap 35l-37j|. 

As T is lowered further, the two sets of spectra evolve 
in a distinctly different fashion. As coherence sets in at 
T c , the spectrum of Fig. QJ, develops the familiar BCS 
coherence peaks at A. As T is reduced below T c the co- 
herence peaks increase in size, while the DoS at E < A 
is suppressed further and seems to vanish as T — > 0. In 
contrast, although a full gap does develop in the spec- 



trum of Fig. [T]d, with a similar magnitude and a van- 
ishing DoS at low energy, the accompanying coherence 
peaks are conspicuously absent. The contrast (below T c ) 
and similarities (above T c ) between the two types of spec- 
tra are highlighted by the corresponding plots of Fig. [TJ; 
andQJi, where more complete T-evolutions of the spectra 
are displayed. 

The similarity between the T-evolution (above T c ), as 
well as the gap magnitude, of both types of spectra in- 
dicate that they share the same physical origin. Because 
the unique shape of the BCS tunneling spectra are known 
to be the result of the Cooper attraction between elec- 
trons, we are led to conclude that the same must be true 
for the locations where the coherence peaks at the gap 
edges are missing from the spectra. 

The appearance of two types of spectra, similar to 
those observed in this work, was reported in the nu- 
merical study of the two-dimensional attractive Hubbard 
model with on-site disorder by Ghosal et al. [3, [l6j |. 
Qualitatively, they found that the regions where coher- 
ence peaks are absent are characterized by rapid fluctu- 
ations of the potential that strongly localize the electron 
states. As a result, they concluded that some of the 
Cooper-pairs that are formed there are unable to partic- 
ipate in the macroscopic coherent state comprising the 
condensate. 

To develop a quantitative understanding of supercon- 
ductivity in this high-disorder limit we turn to more re- 
cent analytic calculations that take into account the frac- 
tal nature of the electronic wave functions close to the 
localization threshold 17, Hij]. In this theory a pair of 



electrons occupying a localized state j experience two 
effects: their net mutual attraction and coherent pair- 
hopping to neighboring states. The observed spectral gap 
A- 7 at a given location contains two distinct contributions 



A^ = Ej{T) 



The first term is a BCS-like exci- 

1/2 



tation energy, determined by E 3 (T) — (ej + hj(T)) 
where Cj is the single electron state energy of the j-th 
state and hj (T) is the local pairing field in this problem. 
In conventional BCS superconductors the pairing field 
is uniform and coincides with the single particle gap. In 
contrast to Ej (T) , the local contribution to the gap A J p is 
not related to the development of a global superconduct- 
ing order-parameter. It results from the Cooper attrac- 
tion between two electrons populating the same localized 
state and it is inversely proportional to the volume of 
that state [3, [13, . Because of the fractality of the 
wave functions the values of A p are large and fluctuate 
strongly for nearly-critical wavefunctions ^j{ r )- An im- 
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mediate conclusion of this theoretical analysis [2! 
that a global superconducting state survives rather deep 
into the localized band, up to the region where an average 
Ap strongly exceeds typical hj(T). In this regime the lo- 
cal pairing field becomes extremely inhomogeneous and 
is characterized by a very broad distribution function. 
The theoretical framework [29|, [38[ explains the T- 
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FIG. 1: Coherent versus incoherent Cooper pairing revealed by local tunneling spectroscopy. T-evolutions of 
the local tunneling tunneling conductance G characterized by the presence (a,c) or absence (b,d) of superconducting coherent 
peaks. The black dashed lines show the spectra measured at T c . a,b, Set of spectra for selected temperatures equal to a fraction 
of the low-T spectral gap. The spectral gap values are A = 560 fj,eV and A = 500 \ieV for a and b respectively (see Methods). 
c,d, Three dimensional view of the same data versus temperature and bias voltage. Spectra of a,c and b,d were measured on 
low and high disorder samples respectively. 



evolutions of the spectra displayed in Fig. [TJ For T > T c , 
global superconductivity is absent (hj = 0) whereas Ap 
remains non-zero indicating the presence of preformed 
Cooper-pairs. This regime is revealed by the pseudogap 
in the DoS shown in Fig. [TJ For T < T c , the develop- 
ment of superconducting correlations at the location of 
Fig. [TJi,c locally gives a non-zero hj , as revealed by the 
progressive growth of the coherence peaks. On the con- 
trary, hj remains nearly zero at the position of Fig.[Tp,d. 
The absence, down to our lowest T, of coherence peaks 
in the gapped DoS at some locations is the fingerprint 
of preformed Cooper pairs that remain localized by the 
strong disorder and do not participate in the condensate. 



Proliferation of localized Cooper pairs when ap- 
proaching the SIT 

We now describe the results of a systematic study of 
our samples. In a disordered system such as ours, a 
quantitative description is provided by the distribution 
functions of the main spectra characteristics obtained at 
different locations on our samples. These characteristics 
are the gap width at low temperatures (defined in meth- 
ods section) and the coherence-peak height defined by 
the ratio R = (G peak - G min )/G min . Here G peak is the 
tunneling conductance at the peak energy and G m i n is 



the minimum tunneling conductance for an energy just 
above the gap. In order to study the dependence on disor- 
der we have collected systematic data from high and low 
disorder samples with T c = 1.2 K and 1.7 K respectively. 
The results, plotted in Fig. O represent 208 I — V traces 
from the high-disorder sample and 2400 traces from the 
low-disorder sample, each trace measured at a different 
location. 

We begin by inspecting the distribution of A. For 
the low disorder sample (Fig. [2Ji), A is normally dis- 
tributed (black solid line in the figure) with an average 
A = 4:90fieV and standard deviation a = 42fieV. The 
increase in disorder (Fig. (2ji) leads to a somewhat lower 
A = 465 /xeV but a larger a — 75/j,eV. This broadening 
of the gap distribution with disorder is similar to that 
observed previously in TiN films [39j . 

A central feature of our results is that the smooth 
evolution of the gap distribution with disorder is ac- 
companied by a very sharp, qualitative, change of the 
coherence-peak height-distribution, V(R), as shown in 
Figs. [2p and[2ji. While the less disordered sample shows 
a well-defined maximum in the V{R) distribution func- 
tion at R 0.3, in the more disordered sample, both 
the maximum of V(R) and most of its weight are shifted 
towards R = 0. This demonstrates that with increasing 
level of disorder, coherent spectra are replaced by those 
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FIG. 2: Disorder-induced inhomogeneity and phase in- 
coherence. Statistics of the local gap value A, (a,c), and the 
local peak-height R = (G pea k — G m i„)/G min , (b,d) (for defi- 
nition see text). Blue (red) histograms correspond to the low 
(high) disorder sample. The difference in the disorder level 
between the two samples is evidenced by the different T c (1.7 
K for low disorder sample and 1.2 K for high disorder sample) 
and normal state resistivity (see Supplementary Fig. SI). 



in which coherent peaks are missing. Close to the SIT, 
this results in a very inhomogeneous state with super- 
conductivity occupying a small fraction of the sample. 
The remaining part of the sample is insulating (although 
the gap in its DoS is also due to Cooper pairing) and its 
relative area grows with increasing disorder. 

Importantly, the fact that incoherently-gapped regions 
proliferate upon the increase of disorder (see Fig. [2]) im- 
plies that the absence of coherence peaks cannot be a 
result of surface contamination as there is no reason to 
expect an abrupt increase of contamination level as the 
SIT is approached. 

Theoretically, the broad statistical distribution of the 
parameters in our samples is a direct consequence of the 
fractality of the wave functions, expected to lead to a 
large and strongly fluctuating values of A p for nearly- 

I. Further, the difference 



critical wavefunctions ipj (r) j 
in the distribution functions displayed in Fig. [2k,, c and 
|2b,d is due to the fact that gap magnitudes are controlled 
by Ap while peak heights are determined by hj . To com- 
pare to the theoretical predictions we show, in Fig.[H the 
computed distribution functions for two levels of disor- 
der. Far from the SIT, the distribution function of hj is 
relatively narrow and it gradually broadens upon increas- 
ing the disorder. In that regime, the typical local field h 
is roughly equal to the average one as shown by the blue 
curve in Fig. [3l Closer to the SIT, the distribution func- 
tion changes dramatically: it acquires a power-law shape 
P(h) ~ h~ a (see red curve in Fig. \S§ with exponent a de- 
creasing upon the growth of disorder (see Supplementary 
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FIG. 3: Local pairing field. Normalized probability distri- 
bution of local pairing amplitudes P(h/T c o) computed in our 
model for slightly different disorders characterized by K = 5 
and K = 6 (see Supplementary Information). The parame- 
ters were chosen so that the transition temperature becomes 
T c0 = W~ 3 E F in the limit of very low disorder and changes 
only by a factor 0.9 between K = 6 and K — 5. We observe 
that in contrast to the moderate effect on T c , the distribution 
functions differ dramatically between these two cases. 



Information). The experimentally observed peak height 
distributions shown in Fig. [2p and d conform well to the 
former and latter situations. 

Despite the inhomogeneous nature of superconductiv- 
ity in our samples, we find that nearly all the spectra are 
characterized by anomalously large 2A/T C (see Fig. |4^). 
The deviation from the classical BCS value can be un- 
derstood by noting that, in tunneling experiments, the 
quantity that is measured is the single-particle DoS and 
the inferred A is the minimal energy for a single electron, 
or hole, excitation. In strongly disordered superconduc- 
tors A is expected to be much larger than the pairing 



field, as observed by Ghosal et al [15|, [16| in their numer 



ical simulations. One therefore expects that, for increas- 
ing disorder, hj will progressively drop to zero, leading 
to a vanishing T c when approaching the SIT while A, 
whose main part is of "incoherent" origin, will remain 
finite and may even increase. In other words, a local 
spectral gap due to Cooper attraction can exist without 
global superconductivity leading in a natural way to the 
large 2A/fcsT c {ks is Boltzmann constant). 



Macroscopic quantum phase coherence probed 
at a local scale 

Finally, to clarify how phase coherence is established in 
such an inhomogeneous superconducting state, we have 
studied the T evolution of several spectra measured at 
different locations r and on samples with different T c . 
For each evolution, we plot in Fig.0^ T pea k(r) /T c , where 
Tpeak (r) is the temperature below which coherence peaks 
start to grow (see Fig. 0p) versus the local 2A(r)/fcsT c . 
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FIG. 4: Onset of the superconducting phase coherence, a, Local onset temperature of the coherence peaks, T pf , a k(r), 
normalized to T c , versus the ratio 2A(r)/fcsT c where A(r) is the low temperature spectral gap. For comparison, we added 
the point T pea k(r) /T c — 1 corresponding to the theoretical BCS ratio 2A/fcgT c = 3.52. b, Thermal evolution of the coherence 
peak height, R (for definition see text), extracted from data of Fig. []J and of the resistivity p of low disorder sample. This plot 
evidences the coincidence between the appearance of the zero-resistance superconducting state at T c with macroscopic phase 
coherence and the onset of the coherence peaks at T pea k. 



While this latter ratio can vary between 6.5 and 11.5, 
Tp e ak(r)/T c remains nearly 1. This proves that peaks 
appear when and only when global superconductivity 
is established regardless of the local gap. In this re- 
spect, our InO films behave very differently from granu- 
lar superconductors where one would expect a constant 
2A(r)//cBT pea fe(r) and a macroscopic T c determined by 
the interplay of the charging and Josephson energies in 
the array of grains 40] . 

The appearance of coherence peaks at the same tem- 
perature where resistance is vanishing (see Fig.0p) allows 
us a direct comparison between local and global measure- 
ments. This is especially striking since both the gap mag- 
nitude and peak heights fluctuate very strongly across 
our samples, while local coherence peaks appear simul- 
taneously at T c . This counter- intuitive behavior allows 
us to rule out a scenario where preformed Cooper-pairs 
would condense locally at different temperatures above 
T c leading to independent superconducting droplets. For 
such an inhomogeneous sample, the zero resistance state 
appearing at T c would simply correspond to the perco- 
lation of these droplets. In a disordered system such as 
ours, because of their fractal nature, the electronic wave- 
functions spread far in the sample and therefore weakly 
overlap with a large number of neighboring states. At T c 
these wavefunctions condense into a single superconduct- 
ing state despite the local fluctuations of the pairing field 
hj. Accordingly, coherence peaks emerge everywhere at 
T c but their height vary from place to place. 

To conclude, we presented the evidence that electrons 
remain Cooper paired at the SIT in disordered InO films. 
The transition is driven by the increase in the number of 
incoherent pairs at the expense of the ones that partic- 
ipate in the condensate. Close to the transition, on the 
superconducting side, A fluctuates strongly indicating a 



spontaneously formed inhomogeneity. The perseverance 
of single electron gap implies that very close to the tran- 
sition on the insulating side, the transport is dominated 
by incoherent Cooper pair hopping. The properties of 
the material in this regime might be very unusual and 
deserve further studies. 



Note that very recent numerical simulations [41( are 
consistent with our interpretation of the data. 



METHODS 

Samples. Our samples are disordered thin films of 
amorphous indium oxide. The films are prepared by us- 
ing electron-beam evaporation of high purity (99.999%) 
I^O^ onto S1O2 in an O2 background. The thicknesses 
of the films studied here are 150 and 300 A as measured 
in situ by a quartz crystal thickness monitor. STM mea- 
surements give a typical rms surface roughness of 1 nm 
on a scanned surface of 1 \xm. Transmission electron mi- 
croscopy studies revealed the amorphous nature of the 
samples without detecting any crystalline inclusion [2lj |. 
In order to perform transport measurements, samples are 
patterned into Hall bridges via a shadow mask and con- 
tacts are made using pressed indium and gold wire. 

Measurements. Transport measurements and tun- 
neling spectroscopy were systematically carried out dur- 
ing the same experiment in a home-built STM cooled 
down to 50mK in an inverted dilution refrigerator. Tem- 
perature of the sample holder, which is weakly coupled 
to the dilution refrigerator, was accurately controlled by 
a R11O2 thermometer and a resistive heater. No mea- 
surable thermal drift of the tip position occurred in our 
experiments between the base temperature of 50mK and 
the highest measured temperature of about 6 K. To per- 
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form spectroscopy, the STM Pt/Ir tip was aligned at the 
center of the sample Hall bridge. The differential con- 
ductance of the tunnel junction, G(V) — was then 
measured by a lock-in amplifier with an alternative volt- 
age modulation of 10 — 30^V added to the ramped bias 
voltage thus allowing to probe the local DoS. The en- 
ergy resolution of the spectroscopy can be described by 
the effective temperature, T e ff — 0.3 K, used to fit the 
superconducting spectra with the theoretical BCS DoS. 
The discrepancy between T e f / and the sample thermome- 
ter is due to unfiltered electromagnetic radiations which 
heat the electrons and from voltage noise generated by 



room temperature electronics 42[. The tunneling cur- 
rent was 0.05 — InA for millivolts bias voltage yielding 
a tunneling resistance in the Mil range, which is much 
higher than the film resistance (< 30 k£l) between the 
STM junction and the contact. This ensures a negligible 
voltage drop across the resistive film in series with the 
STM junction during spectroscopy above T c . For statis- 
tics analysis of the gap distribution we used the practical 
definition of the gap A/e = l.lV ma x, where e is the elec- 
tron charge and V max was defined as the bias voltage at 
which the numerically computed dG(V)/dV is maximum 
at the gap edge. The factor 1.1 compensates the shift of 
Vmax to lower values that is due to thermal broaden- 
ing of the differential tunneling conductance. This rela- 
tion which is valid for superconducting spectra was also 
used to quantify the gap value for spectra without peaks. 
Transport measurements of the films (see Supplementary 
Information) were carried out with a low frequency lock- 
in amplifier technique in a four terminal configuration 
with excitations currents of InA below 5K and of lOnA 
at higher temperatures. 



* Present address: Institut Neel, CNRS and Universite 
Joseph Fourier, BP 166, 38042 Grenoble, France 
[1] Goldman, A. M., and Markovic, N. Superconductor- 
insulator transitions in the two-dimensional limit. Phys. 
Today. 51, 39 (1998). 
[2] Bezryadin, A., Lau, C. N., Tinkham, M. Quantum sup- 
pression of superconductivity in ultrathin nanowires. Na- 
ture 404, 971-974 (2000). 
[3] Steiner, M. A., Boebinger, G., Kapitulnik A. Possible 
Field- Tuned Superconductor-Insulator Transition in High- 
Tc Superconductors: Implications for Pairing at High 
Magnetic Fields. Phys. Rev. Lett. 94, 107008 (2005). 
[4] Sanchez-Palencia, L., Lewenstein, M. Disordered quantum 
gases under control. Nature Physics 6, 87-95 (2010). 
Basko D. M., Aleiner, I. L., Altshuler, B. L. Metal- 
insulator transition in a weakly interacting many-electron 
system with localized single-particle states. Ann. Phys. 
321, 1126-1205 (2006). 
[6] Oganesyan, V., Huse, D. A. Localization of interacting 
fermions at high temperature. Phys. Rev. B 75, 155111 
(2007). 

[7] Gornyi, I. V., Mirlin, A. D., Polyakov, D. G. Interacting 



[5 



Electrons in Disordered Wires: Anderson Localization and 
Low-T Transport. Phys. Rev. Lett. 95, 206603 (2005). 

[8] Anderson, P.W. Theory of dirty superconductors. J. Phys. 
Chem. Solids 11, 26-30 (1959). 

[9] Abrikosov, A. A., Gorkov, L. P. On the theory of supercon- 
ducting alloys .1. The electrodynamics of alloys at absolute 
zero. Sov. Phys. JETP 8, 1090 (1958). 

[10] Bardeen, J., Cooper, L. N., Schrieffer, J. R. Theory of 
Superconductivity. Phys. Rev. 108, 1175-1204 (1957). 

[11] Bulaevskii, L. N., Sadovskii, M. V. Localization and su- 
perconductivity. Pisma ZhETF 39, 524 (1984) [Sov. Phys. - 
JETP Letters, 39, 640-643 (1984) ]. Sadovskii, M. V. Su- 
perconductivity and localization. Phys. Rep. 282, 226-348 
(1997). 

[12] Kapitulnik, A., Kotliar, G. Anderson Localization and 

the Theory of Dirty Superconductors. Phys. Rev. Lett. 54, 

473-476 (1985). 
[13] Ma, M., Lee, P. A. Localized superconductors. Phys. Rev. 

B 32, 5658-5667 (1985). 
[14] Kotliar, G., Kapitulnik, A. Anderson localization and the 

theory of dirty superconductors. II Phys. Rev. B 33, 3146- 

3157 (1986). 

[15] Ghosal, A., Randeria, M., Trivedi, N. Role of Spatial 
Amplitude Fluctuations in Highly Disordered s-Wave Su- 
perconductors. Phys. Rev. Lett. 81, 3940-3943 (1998). 

[16] Ghosal, A., Randeria, M., Trivedi, N. Inhomogeneous 
pairing in highly disordered s-wave superconductors. 
Phys. Rev. B 65, 014501 (2001). 

[17] Feigel'man, M. V., Ioffe, L. B., Kravtsov, V. E., 
Yuzbashyan, E. A. Eigenfunction Fractality and Pseudo- 
gap State near the Superconductor-Insulator Transition. 
Phys. Rev. Lett. 98, 027001 (2007). 

[18] Finkel'stein, A. M. Superconducting transition tempera- 
ture in amorphous films. Pis'ma Zh. Eksp. Teor. Fiz. 45, 
37 (1987) [Sov. Phys. JETP Lett. 45, 46 (1987)]. 

[19] Haviland, D. B., Liu, Y., Goldman, A. M. Onset of su- 
perconductivity in the two-dimensional limit. Phys. Rev. 
Lett. 62, 2180-2183 (1989). 

[20] Hebard, A. F., Paalanen, M. A. Magnetic-field-tuned 
superconductor-insulator transition in two-dimensional 
films. Phys. Rev. Lett. 65, 927-930 (1990). 

[21] Shahar, D., Ovadyahu, Z. Superconductivity near the 
mobility edge. Phys. Rev. B 46, 10917-10922 (1992). 

[22] Gantmakher, V. F., Golubkov, M. V., Dolgopolov, V. 
T., Tsydynzhapov, G. E., Shahskin, A. A. Destruction 
of localized electron pairs above the magnetic-field-driven 
superconductor-insulator transition in amorphous In-0 
films. JETP Letters 68, 363-369 (1998). 

[23] Steiner M., Kapitulnik, A. Superconductivity in the 
insulating phase above the field-tuned superconductor- 
insulator transition in disordered indium oxide films. Phys- 
ica C 422, 16-26 (2005). 

[24] Crane, R. et al. Survival of superconducting correlations 
across the two-dimensional superconductor-insulator tran- 
sition: A finite-frequency study. Phys. Rev. B 75, 184530 
(2007). 

[25] Hadacek, N., Sanquer, M., Villegier, J.-C. Double reen- 
trant superconductor-insulator transition in thin TiN 
films. Phys. Rev. B 69, 024505 (2004). 

[26] Baturina, T. I., Mironov, A. Yu., Vinokur, V. M., Bak- 
lanov, M. R., Strunk, C. Localized Superconductivity 
in the Quantum-Critical Region of the Disorder-Driven 
Superconductor-Insulator Transition in TiN Thin Films. 
Phys. Rev. Lett. 99, 257003 (2007). 



7 



[27] For yet a third approach, see Dubi, Y., Meir, Y., Avishai, 
Y. Nature of the superconductor-insulator transition in 
disordered superconductors. Nature 449, 876 (2007). 

[28] Fisher, M. P. A. Quantum phase transitions in disor- 
dered two-dimensional superconductors. Phys. Rev. Lett. 
65, 923-926 (1990). 

[29] Feigel'man, M. V., Ioffe, L. B., Kravtsov, V. E., Cuevas, 
E. Fractal superconductivity near localization threshold. 
Annals of Physics 325, 1390-1478 (2010). 

[30] Sambandamurthy, G., Engel, L. W., Johansson, A., Sha- 
har, D. Superconductivity- Related Insulating Behavior. 
Phys. Rev. Lett. 92, 107005 (2004). 

[31] Baturina, T. I., Strunk, C, Baklanov, M. R., Satta, 
A. Quantum Metallicity on the High-Field Side of the 
Superconductor-Insulator Transition. Phys. Rev. Lett. 98, 
127003 (2007). 

[32] Stewart, Jr., M. D., Yin, A., Xu, J. M., Valles, Jr., J. 
M. Superconducting Pair Correlations in an Amorphous 
Insulating Nanohoneycomb Film. Science 318, 1273-1275 
(2007). 

[33] Nguyen, H. Q. et al. Observation of Giant Positive Mag- 
netoresistance in a Cooper Pair Insulator. Phys. Rev. Lett. 
103, 157001 (2009). 

[34] Giaever, I. Energy Gap in Superconductors Measured by 
Electron Tunneling. Phys. Rev. Lett. 5, 147-148 (1960). 

[35] Fischer, 0., Kugler, M., Maggio-Aprile, I., Berthod, 
C, Renner, C. Scanning tunneling spectroscopy of high- 
temperature superconductors. Rev. Mod. Phys. 79, 353- 
419 (2007). 

[36] Timusk, T., Statt, B. The pseudogap in high- 
temperature superconductors: an experimental survey. 
Rep. Prog. Phys. 62, 61-122 (1999). 

[37] Sacepe, B. et al. Pseudogap in a thin film of a conven- 
tional superconductor. Nat. Commun. 1, 140 (2010). 

[38] Ioffe, L. B., Mezard, M. Disorder-Driven Quantum Phase 
Transitions in Superconductors and Magnets. Phys. Rev. 
Lett. 105, 037001 (2010). 

[39] Sacepe, B. et al. Disorder-Induced Inhomogeneities of 
the Superconducting State Close to the Superconductor- 
Insulator Transition. Phys. Rev. Lett. 101, 157006 (2008). 

[40] Efetov, K. B. Phase-transition in granulated supercon- 
ductors. Zh. Eksp. Teor. Fiz. 78, 2017-2032 (1980) [Sov. 
Phys. JETP 51, 1015 (1980)]. 

[41] Bouadim, K., Loh, Y. L., Randeria, M. Trivedi, N. Single 
and two-particle energy gaps across the disorder-driven 
superconductor-insulator transition. arXiv:1011.3275 

[42] Le Sueur, H., Joyez, P. Room-temperature tunnel cur- 
rent amplifier and experimental setup for high resolution 
electronic spectroscopy in millikelvin scanning tunneling 
microscope experiments. Rev. Sci. Instrum. 77, 123701 
(2006). 



Acknowledgements We thank A. Finkel'stein, V. 
Kravtsov, M. Mezard, Z. Ovadyahu and N. Trivedi for 
valuable discussions. D.S. and M.O. acknowledge the Is- 
raeli Science Foundation and the Minerva Fund. 



Author Contributions B.S. and T.D. equally con- 
tributed to this work. 



SUPPLEMENTARY INFORMATION 
Transport measurements 

In this section we present the transport properties of 
our highly disordered amorphous InO films. The T- 
cvolutions of the resistivity for the three samples studied 
in this work are shown in Fig. IS5I All films present an 
insulating trend: The resistivity increases as T is low- 
ered from room temperature. Despite this clear signa- 
ture of electron localization, the films undergo at low 
temperature a superconducting transition with critical 
temperature 1 — 2K (see inset of Fig. IS5[) . The con- 
tinuous decrease of the critical temperature upon in- 
creasing disorder is the hallmark of the disorder-driven 
superconductor-insulator transition that occurs in our 
InO films. 




low disorder 
intermediate disorder 
high disorder 



50 100 150 200 250 300 
T[K] 



FIG. S5: Resistivity versus temperature of amorphous 
indium oxyde films. Thicknesses of the samples are 150 
A for low disorder sample and 300 A for high and interme- 
diate disorder samples. Inset: superconducting transitions 
measured during the cooling of the STM setup. The critical 
temperatures defined as the vanishing resistance (see main 
text) are 1.7 K, 1AK and 1.2 K for low, intermediate and 
high disorder sample respectively. 



Summary of the theory of superconductivity 
close to the mobility edge 

Here we sketch the main ingredient of the theory de- 
veloped in 0-[3l Q and its application to the data. In 
the range of strongly developed pseudogap Ap » T c 
electron "orbitals" -0j ( r ) are populated by either zero or 
two electrons at low T's. In this case the whole Hilbert 
space of the electronic problem reduces to the "pseu- 
dospin" subspace described by operators formally equiv- 
alent to spin-i variables S. , 5 J associated with each or- 
ci~ ,<& ■ is the pair creation operator and 



bital: S+ 

2Sj = jC^-j + c^.cu — 1; this representation was pro- 
posed originally by Anderson in 4]. The development 
of superconducting coherence is described by the pseu- 
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dospin Hamiltonian 

H PS = 2 ]T - f £ M« (S+S, + Sr S +), (1) 

j ij 

with matrix elements My = J d 3 r0 2 (r)-0|(r) and cou- 
pling constant On-site energies Cj are distributed over 
wide band with bare DoS isq. The pseudogaped regime 
Ap > T c is realized [l|, Q when typical level spacing 
inside localization volume, 5 ~ l/i>o£f oc , i s the l ar gest 
energy scale, S 3> Ap. The local pairing fields /ij(T) de- 
termine the average transverse pseudospin component via 
(Sf) = (hi/2E i )tauh(E i /T); they obey self-consistency 
equation 

fti = §E M «F tM,h §' E i = J h S +f %' (2) 

At very large 5, long-range coherence is not estab- 
lished even at T = 0, and insulating ground state takes 
over. The quantum phase transition between the ordered 
state with non-zero (S x ' v ) ^ and disordered one with 
{S x,v ) — was studied in Ref. 0, Q within the simpli- 
fied model defined by the Hamiltonian ((T|) on the Bethe 
lattice with coordination number K and all nonzero cou- 
plings Mij = 1/ K. The model is characterized by dimen- 
sionless coupling A = gv$ <C 1. The result of this study 
is that in the limit K > K x = Ae 1 / 2A usual BCS-like 
mean-field theory works well and uniform superconduct- 
ing state occurs below T c o = is ~ 1 e~ 1 / x . In particular, 
below T c q nonzero thermal average hj(T) appear. 

Fluctuations become important at K < K\ = Ae 1 ^ 2A . 
They lead to two effects: critical temperature T C (K) start 
to drop down with K decrease and the order parameter 
hi (T) becomes strongly inhomogeneous as function of the 
site number i and cannot be expressed in terms of local 
energy ej only. Eventually, transition temperature van- 
ishes at K = K2 = Ae 1 / eA , the ground-state at K < K2 
is a many-body insulator, with discrete spectrum of low- 
lying excitations. 

In the whole region K2 < K < K\ on the supercon- 
ducting side of quantum phase transition, nonzero order 
parameter hj (T) appears below well-defined global tran- 
sition temperature T C {K) and leads to the growth of co- 
herence peaks seen in the tunneling spectra, Fig. la,c. At 
the same time, the relative heights of these peaks are pro- 
portional to the values of hj(T) for the electron orbitals 
which have considerable weights near the tip posi- 
tion r. The values of hj fluctuate strongly between differ- 
ent sites j with very close local energies €j . The strength 



of these fluctuations is characterized by the distribution 
function P(h) with a long tail: close to the transition 
line T C (K) the distribution is P(h) - (l/h ){h Q /h) a in 
a wide range of h/hg. The exponent a decreases with 
decrease of K, i.e. in physical terms, with increase of 
disorder. In particular, at K2 < K < K\ the exponent 
1 < a < 2 was found [3, 5], indicating that simple av- 
erage h av = (h) and typical value h typ — exp{ (log (h))) 
differ qualitatively, with h typ <C h av ( h typ character- 
izes behaviour of almost any specific sample, whereas h av 
corresponds to a contribution of extremely rare fluctua- 
tions in the ensemble average over many samples). We 
show in Fig. 3 theoretical results for two nearby values 
of the coordination number K — 5 and K = 6, both at 
the dimensionless coupling constant A = 0.128. While 
the values of the critical temperature are close in these 
two cases (T C (K = b)/T c {K = 6) « 0.9), the shape of 
P(h) distribution changes considerably: the "tail" with 
h <C htyp becomes much more pronounced for the K = 5 
case that corresponds to stronger disoder. This theoret- 
ical result is in good agreement with our observations 
shown in Fig. 2b, d: distribution of peak heights is mod- 
erately narrow in the less disordered sample shown in 
Fig. 2b and becomes very broad (like the P(h) distribu- 
tion function at K < K\) for the more disordered sample 
in Fig. 2d. 
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